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ABSTRACT 

We analyze Spitzer GLIMPSE, MSX, and WMAP images of the Milky Way 
to identify 8 micron and free-free sources in the Galaxy. Seventy-two of the 
eighty-eight WMAP sources have coverage in the GLIMPSE and MSX surveys 
suitable for identifying massive star forming complexes (SFC). We measure the 
ionizing luminosity functions of the SFCs and study their role in the turbulent 
motion of the Galaxy's molecular gas. We find a total Galactic free-free flux fi, 
= 46177.6 Jy; the 72 WMAP sources with full 8 micron coverage account for 
34263.5 Jy (~ 75%), with both measurements made at z/=94GHz (W band). 
We find a total of 280 SFCs, of which 168 have unique kinematic distances and 
free-free luminosities. We use a simple model for the radial distribution of star 
formation to estimate the free-free and ionizing luminosity for the sources lacking 
distance determinations. The total dust-corrected ionizing luminosity is Q = 
2.9 ± 0.5 X 10^^ photons s^^, which implies a galactic star formation rate of 
= 1.2 ± 0.2 Mq yj:-\ We present the (ionizing) luminosity function of 
the SFCs, and show that 24 sources emit half the ionizing luminosity of the 
Galaxy. The SFCs appear as bubbles in GLIMPSE or MSX images; the radial 
velocities associated with the bubble walls allow us to infer the expansion velocity 
of the bubbles. We calculate the kinetic luminosity of the bubble expansion and 
compare it to the turbulent luminosity of the inner molecular disk. SFCs emitting 
80% of the total galactic free-free luminosity produce a kinetic luminosity equal 
to 65% of the turbulent luminosity in the inner molecular disk. This suggests 
that the expansion of the bubbles is a major driver of the turbulent motion of 
the inner Milky Way molecular gas. 
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Introduction 



Turbulence is ubiquitous in astrophysical fluids, from the intergalactic medium of galaxy 
clusters through the interstellar medium of galaxies to protostellar accretion disks. Obser- 
vations of CO linewidths corresponding to supersonic velocity dispersions have led to the 
suggestion that turbulence in molecular clouds m ay slow gravitational contraction and limit 
the rate of star formation, e.g., IShu et al.l ( 119871 ). If the motions are turbulence-dominated, 
there must be a driving energy source, since undriven turbul ence decays on the order of a 
crossing time (IGoldreich fc Kwanlll974j : Mac Low et al.lll998l ). Currently proposed sources 
can be divided into tw o groups: gravitational energ y released by accretion onto or through 
the galactic disk (e.g., iKlessen fc Hennebelld 120101 ). and nuclear energy. Gr avitational en- 



ergy can generate turbulence via instabil ities such as bars or spiral arms 



on large scales, and protostellar outflows (ILi fc Nakamurall2006l : IWang et al. 



Li et all 120051 ) 
2010[ l on small 



scales, while nuclear e nergy couples t o the interstellar me dium (ISM) turbule nce via ex- 
pandin g H II regions (IMatznerl 120021) . radiation pressure (IMurray et al.l 120051 ) . or super- 
novae (IMac Low fc KlessenI I2OO4J ) . In the absence of star formation activities such as in the 
extended H I disk, magnetorotational instab ility is proposed as one of the sources of tur- 



bulent motion (e.g. ISellwood fc Balbuslll999l ). None of the proposed mechanisms, however. 



have strong observational support. 



Recently. [Murray fc Rahman! (120101 . hereafter MRIO) identified 88 Wilkinson Microwave 
Anisotropy Probe (WMAP) galactic free-free sources in an effort to determine the galactic 
star formation rate (SFR). These WMAP sources contain bubble-like structures with radii 
ranging from 5 to 100 pc. Since the 18 most luminous sources identified produced over half 
the total ionizing luminosity of the Galaxy, they proposed that these sources contain over 
half the Galactic O star population, most likely in massive cluster^ The ionizing photons 
produced by the clusters leak out of their host bubbles and ionize the surrounding area. 
Th ese leak ii ig ph otons create the extended low density (ELD) regions, originally described 
by iMezgerl (119781 ) as the dominant source of galactic radio free-free emission; the WMAP 
sources are argued to be the ELDs. 

Noting the existen ce of numerous known H II regions on the shells or walls of the bubbles 
in the WMAP sources, iRahman fc Murray! (!2010! . hereafter RMIO) identified 40 star forming 
complexes (SFCs) in the 13 most luminous WMAP sources, responsible for one-third of the 
star formation in the Milky Way. The kinetic luminosity due to the expansion of the bubbles 
associated with the SFCs can account for a third of the turbulent luminosity in the Milky 
Way's molecular gas disk, suggesting that the bubble expansion is the primary driver of 
turbulence. They also found that the initial driver of the bubbles was not supernovae, as the 
regions were selected by their free-free luminosity; the free-free luminosity of a star cluster 



''A candidate OB association has been found in G298 a nd a follow-up spectro scopic analysis confirms 15 
0-type and 3 luminous blue variables or Wolf-Rayet stars ( Rahman et al. 2011a! !b). 
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fades rapidly after the first supernovae explode, because the most massive stars dominate 
the free-free emission. Since the SFCs selected by RMIO are in the most luminous WMAP 
sources, further investigation is required to determine if the relationship between free-free 
and kinetic luminosity scales to the less luminous sources. 

In this paper, we perform a follow-up study by investigating the remaining 75 WMAP 
sources and revisiting the previously-studied 13. In ^ we describe the selection of the 
SFCs, discuss their expansion velocities, distance and free-free flux determination, and show 
comparison statistics to RMIO. In ^ we refine the galactic ionizing luminosity and star 
formation rate, present an empirical relationship between the 8/im and free-free emission, 
and show that the expansion of the SFCs is a significant driving mechanism of turbulent 
motion of the molecular gas inside the solar circle. We discuss the results in §11 In §3 we 
summarize the results. 



2. Data Analysis 

2.1. Selection of Star Forming Complexes 

MRIO identified 88 galactic free-free sources in the WMAP maximum entropy method 
free-free foreground emission map: 87 in the Galaxy, the 88th being the LMC. The 13 most 
luminous WMAP sources were then analyzed by RMIO using radial velocity measurements 
taken from the literature, combined with the morphological features in 8/im images. Here, 
we investigate all the WMAP free-free sources to obtain a complete sample of the Galactic 
free-free emitting SFCs. 

In order to encompass an entire WMAP source in a single image, band 4 (8 /im) mosaics 
from the Spitzer Galactic Legacy Infrared Midplane Survey Extraordinaire (GLIMPSE) are 
stitched together using the Montage package. Forty-nine sources do not have sufficient 
GLIMPSE coverage; for these sources, Midcourse Space Experiment (MSX) band A images, 
also at 8 /xm but with lower sensitivity and spatial resolution, are used instead. Visual 
inspection of the 8/im images revealed bubble-like morphologies outlined by bright regions 
along the walls that are well correlated with H II regions found from SIMBAD. Catalogues 
of radio recombination line (RRL) measurements taken from the literature are used to find 
the radial velocities of these H II regions. In cases where there are no RRL measurement. 
Ha or molecular line velocities are used (non-RRL velocities are used only for SFC selection 
and are excluded in further analysis). SFCs are selected based on the visual inspection of 
the shell morphology in the 8/im images, and a small range in RRL velocities, where we set 
th e maximum spread to be 15km s~^, based on the theoretical study of expanding bubbles 



by iHarper-Clark fc Murrayl (120091 ). Thus, we interpret the dispersion in the velocities of the 



H II regions as the expansive motions of the bubbles inside the SFC. 
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Often 2 or 3 additional SFCs are found just outside the WMAP ellipses from MRIO; in 
these cases, the WMAP source is either resized and reoriented to encompass all the observed 
SFCs, or new WMAP sources are created if RRL velocities of SFCs outside the original 
WMAP ellipses are significantly different from those inside. On the other hand, some SFCs 
are found to lie across multiple WMAP sources; in these cases, WMAP sources are merged 
appropriately. The mergers effectively eliminate three MRIO sources while the splits add 4 
new sources. Both the original and the new WMAP sources are shown in Figure [H 

MRIO found 87 free-free sources in the Galaxy. After merging and splitting sources, 
we have 88 in which we can potentially find star forming complexes. However, ten of our 
WMAP sources do not have any 8/im coverage and are excluded from our SFC search. After 
the further exclusion of the WMAP regions with only partial 8/im coverage (G64, G107, 
G118) or a lack of bubble detection (G87, G130, G271), a total of 72 WMAP sources are 
suitable for finding SFCs. The total free-free flux of the 72 WMAP sources is 34263.5 Jy, 
which is 74% of the total Galactic free-free flux 46177.6 Jy. 

We find a total of 280 SFCs in these 72 sources; a sample selection is shown in Figure 

m 

We note that MRIO estimate the total galactic free-free flux as 54211.6 Jy, measured 
in Ka band (33GHz). By summing up all the pixels in each of the 5 WMAP band images 
and plotting as a function of the frequency, we find Fy oc ^-o-i6± o.ii^ using a 12% error as 



estimated from the Ha emission discrepancy (IBennett et al.ll2003l ). This is consistent with 



the frequency dependence of free-free emission, due primarily to the variation in the Gaunt 
factor (see equation [1]) . It is this weak frequency dependence which leads to the lower free- 
free flux quoted in this paper. All the free-free fluxes quoted in this paper are measured in 
the W band at 94 GHz from the WMAP free- free foreground emission map (MRIO). 



2.2. Properties of Star Forming Complexes 

2.2.1. Radial Velocities and Distances 

We interpret the morphology seen in the 8/im images as the result of quasi-spherical 
expansion, assuming that the RRL velocities sample the expansion rate of the bubble wall. 
However, due to the line-of-sight projection of the three-dimensional expansion motion and 
limited sampling of RRL velocities within the SFCs, a geometric correction is required to 
convert the measured velocity spread into an expansion velocity. The method outlined in 
RMIO is used, in which the measured half-spread velocity (defined as {vmax — Vmin)/1) is 
multiplied by a geometric factor found by running a Monte Carlo simulation where each 
velocity measurement is de-projected by assigning a random pair of spherical coordinate 
angles. The mean correction factor is 2.0, 1.4 and 1.2 for 5, 10 and 15 velocity measurements 
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Fig. 1. — Histogram-equalized image of all pixels above 0.02Jy in W band at 94 GHz from 
the WMAP free-free foreground emission map with the WMAP sources overlaid (MRIO). 
Original positions and sizes are indicated by red ellipses and the resized versions are indicated 
by turquoise ellipses. 




Fig. 2.— Spitzer GLIMPSE image of a part of the G317 WMAP source with SFCs indicated 
by blue elhpses and the H II region velocities indicated by red circles. 
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respectively (RMIO). 



Defining the median of tlie local standard of rest (LSR) velocities of the H II (RRL) re- 
gions as the velocity of the host SFC, a kinematic distance is determined by fitting th e galactic 
l ongit ude of the centroid of the SFC and the LSR velocity of the SFC using the I Clemens 
( 119851 ) rotation curve, assuming that the distance from the sun to the Galactic center is Rq 
= 8.5 kpc, and a circular velocity 9 = 220 km s~^. The kinematic distance of SFCs with 
highly peculiar (forbidden) velocities are re-calculated allowing for a deviation in the velocity 
of ± 20 km s-\ 

SFCs that lie in the inner galaxy (r < Rq) suffer from the kinematic distance ambiguity 
(KDA); the source could lie at either the near or far points where the projected line of sight 
velocity matches the measured velocity. Absorption line data from the literature are used 
to break the degeneracy. We use absorption line data from a number of sources. First, 
Anderson &: Banial (120091 ) provide both H I emission/absorption and self-absorption data, 
and a corresponding estimate of the reliability of the result. In case these two measurements 
do not agree, we use the result with the higher assigned confidence level. In case the methods 
disagree and are given the same level of confidence, we compare their quoted distance to our 
distance derived from RRL velocities. 



Other sources from the radial velocity litera ture that we e mploy include iBusfield et al. 



(120061), wh o use H I s elf-ab sorpt ion data whi^ 



data. Both lFish et al.l ( 120031 ) and lKolpak et al 



Sewilo et al.l ( 120041 ) use H2CO absorption 



(l2003h _ jrovi de H I line absorption (measured 



against H II continuum background). ICaswell fc Hayned (119871 ) use a combination of H I 



and H2CO absorption data as well as optical counterparts to the H II regions. 



Palagi et al.l (Il993[ ) and I Gillespie et al.l (Il977l ) quote only the distances (and not the 



radial velocity fr), where the former re s olve t he KDA by comparing with the result of 
Downes et al.l (Il980[ ) and ISolomon et al.l (Il987| ) and the latter by using H2CO absorption 
data and the spectrophotometry of associated stars. In these and similar cases, when only 
the absorption line-based distance is given (and not Vr), we compare the distance in the 
literature to the distance determined from the radio recombination line Vr combined with 
the Clemens rotation curve, and choose the closest match. 

Given a SFC with multiple resolutions of the KDA, we assume majority rules; for 
example, if a SFC has 5 KDA measurements with 3 near and 2 far, we take the near distance. 
In case of a tie, we take the measurement with (literature reported) higher confidence. Failing 
that, we check the vertical distance of the SFC from the plane and take the near distance if 
the far distance gives the height greater than lOOpc. If the ambiguity remains, we discard 
measurements from regions that do not project onto the shell walls in 8/im images. Finally, 
we check the physical size of the SFC — based on the hypothesis that the bubbles are blown 
by massive stars, SFC size can be no bigger than the expansion velocity times the lifetime of 
massive O stars (4 Myr). There are 26 SFCs with Tdyn = physical radius/ expansion velocity 
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> 4 Myr. Of these 26 SFCs, 7 of them do not have corrected mean-spread velocity. As the 
median discrepancy of the dynamical time is ~ 49%, not a factor of 2 or higher, and the 
physical radius of the bubble is on the order of or less than 100 pc, we interpret them as an 
underestimate of expansion velocity. 

Three SFCs are assigned a kinematic distance manually: G43 SFC92, G359 SFCO, G359 



SFCl . G43 SFC92 is assigned 6 kpc from the tangent point measurement of lAnderson fc Bania 



( I2OO9I ) (the measured near and far distance are 5 kpc and 7 kpc respectively). G359 SFCO 
and SFCl are both assumed to lie at 8 kpc as they are most likely associated with Arches 
cluster and Quintuplet cluster, respectively. 

Use of the Galactic rotation curve as a distance estimator fails near the Galactic centre, 
giving erroneously far or near distance for SFCs with small Galactic longitude. As will be 
described in ^2.2.41 we caii infer the association of SFCs with the ~ 3 kpc molecular ring 
(see, eg., iDame et al.l 1200 ll ) using only Galactic longitude and LSR velocity. If a SFC we 



expect to be part of the molecular ring is assigned a distance larger than 14 kpc , we manually 
change the distance to 12 kpc ; conversely, if the SFC is assigned to a distance closer than 
3 kpc , we manually change the distance to 4 kpc . 

A total of 194 SFCs require KDA resolution; we find unique distances for 113 of these 
194 sources. Combined with the SFCs that do not require the resolution, we identify a 
unique distance estimate to 168 out of our 280 SFCs. 



Reid et al.l (l2009h . based o n their parallax s tudy, suggest that 6 = 250 km s ^ rather 
than 220 km s~^. Assuming that lReid et al.l ( 120091 ) is correct, using the distance calculated by 



the IClemensI ( 1l985l ) curve results in outer galaxy distances being, on average, overestimated 
by a factor of 1.5. The inner g alaxy source distance s aremuch less sensitive to 6', changing by 
approximately 1% on average. iMcMillan fc Binneyl ( I2OIOI ) point out that Rq and 9 cannot be 
reliably determined an d are highly sensi tive to model parameters; furthermore, they point 
out that the result of iReid et al.l (120091 ) may not be accurate due to the underestimated 
peculiar motion of the Sun with respect to the LSR. Nevertheless, our results shows that the 
distance to the inner galaxy SFCs is largely insensitive to the value of 9, and the distance 
discrepancy in the outer sources, as will be described in § 13. 2[ can alter our SFR by ~ 6%. 



2.2.2. Free-free fluxes 

The WMAP beam is significantly larger than the typical SFC radius, resulting in spatial 
confusion between the free-free emission of different SFCs. To overcome this confusion and 
attribute free-free flux estimates to individual SFCs, we divide the total WMAP flux among 
the corresponding SFCs by the relative 8/im flux. (Here we define the total WMAP free-free 
flux as the isophotal flux calculated in MRIO). The relationship between free- free flux and 
8/im emission is discussed in §3.31 where it is shown that 8/im is a good tracer of free-free 
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emission. Previously, iMurrayi (1201 ih used relative RRL flux to divide up total WMAP free- 
free flux. These two methods agree well, with correlation coefficient ~0.7. We prefer using 
the 8/im flux because it does not require uniform sampling (both in number and wavelength) 
of RRL measurements among the SFCs in the same host WMAP, sampling which is rarely 
available. All 280 SFCs are assigned a free-free flux using this method. As a result, we find 
the total flux of the SFCs to be 34234.8 Jy (the 0.1% difference from the total 72 WMAP 
source free-free flux appears to be due to rounding error in our python scripts), leaving 26% 
of the total Galactic free-free flux 46177.6 Jy unaccounted for. Ultimately, 166 SFCs have 
both unique kinematic distance and free-free flux determinations. 



2.2.3. Radial velocity and distance uncertainties due to streaming motions of H II gas 

It is possible that the H II region line-of-sight velocities may probe the superposition of 
the streaming motion of the ionized ga s driven locally by the O stars powering the free-free 



emission (~ 10 km s~ iMatzned 120021 ) . in combination with the expansion of the bubble. 
However, Tabled] shows that the average dynamical time of the SFCs Tdyn = R/Avc is 4 Myr. 
This is roughly the lifetime of the most massive O stars, which dominate the production 
of ionizing radiation. In other words, the bubble expansion velocity, calculated from the 
bubble size divided by the upper limit for the age of the system is of order 10 km s~^. This 
is consistent with the independent measurement of the expansion velocity provided by the 
RRL data: the mean expansion velocity of all the SFCs is found to be 12 ± 8 km s~^. 

It is also possible that the calculated LSR velocity of the SFC may be in error, affecting 
the kinematic distance measurement. Figure |3] shows the error in distance AD as a function 
of kinematic distance D. The error contribution due to the (assumed) streaming motion of 
ionized gas is similar to the error due to the velocity spread inferred from (multiple) radio 
recombination line velocities. While there are a few cases where the fractional error is ~ 10 
or even higher, the majority of the SFCs have a 36% error in the kinematic distance. 

As a sanity check, the LSR velocity and t he expansion velocity of SFC 264 in G353 can be 



compared to the results of ICappa et al.l (l201l[ ). who find an expansion velocity ~ 5— 9 km s~^ 
which is in approximate agreement with the result quoted in this paper, 6.7 km s~^. Using 
a photometric distance to the cluster Pismis 24, which appears to power the HII region, they 
find D = 2.5 kpc, substantially larger than our D = 0.9 kpc. Photometric distances suffer 
from uncertainties in reddening, while kinematic distance estimates suffer from uncertainties 
in the velocity measurement and the accuracy of the rotation curve. For example, the 
two distance estimates could be reconciled by correcting the LSR velocity —3.8 km s~^ for 
a ~ 10 km s~^ streaming motion of the ionized gas, yielding Vr = —13.8 km s~^ and a 
kinematic distance measurement of 2.9 kpc. We conclude, therefore, that the two distance 
results agree within the uncertainties. We note that a study of CO map or any neutral gas 
tracer will provide a better estimate of both the expansion velocity and central velocity of 
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the bubbles. 



2.2.4- Uncertainties due to Spiral Arms 



While the aforementioned derivation of kinematic distance assumes circular motion of 
the SFCs, the kinematics of Milky Way interstellar medium (ISM) are af f ected by the non- 
circular motions of spiral structures (IBinney fc Merrifieldlll998l ). iGomezI (120061 ) shows that 
the kinematic distances derived from a measured rotation curve have errors from 0.5 kpc to 
3kpc near the spiral arms. Takin g 1.5 kpc as t he typical error in the kinematic distance 
near the Galactic spiral arms (see iGomezI |2006| . their Figure 5), we add this error to the 
kinematic distance error propagation of SFCs that are likely associated with spiral arms. 
Figure H] shows the line of sight ra dial velocity versus the Galactic longitude of the SFCs. 
Comparing Figure H] to Figure 3 of iDame et al.l ( l200ll ). we associate any SFCs that satisfy 
any of the following criteria with a spiral arm: 



Local Arm: 100° ^ / ^ 160° and \vlsr\ ^ 20 km s~^ 
Perseus & Outer Arm: 30° < / ^ 180° and Vlsr < -20 km s"^ 
Carina & Perseus Arm: 180° ^ / ^ 330° and vlsr > 20 km s"^ 
Molecular Ring: |/| ^ 30° and \vLSR,rot\ = 20km s^^ 



where VLSR,rot corresponds to vlsr rotated in the y-1 plane such that the molecular ring 
structure in Figure H] is aligned along the vlsr = km s~^ axis. We include the molecular 
ring in the uncertainty analysis as w ell since it is likel y composed of multiple spiral arms 
rooted from the Galactic central bar (IDame et al.l 120011 ). 



2.2.5. Catalogue of SFCs 

The complete SFC catalogue is presented in Tabled] with the columns as follows: Column 
1 gives the catalogue number, Column 2 the host WMAP source number, Columns 3 and 4 
galactic longitude and latitude. Columns 5 and 6 semi-major and semi-minor axes. Column 
7 position angle. Column 8 mean LSR velocity. Column 9 measured half-spread velocity. 
Column 10 estimated radial expansion velocity. Column 11 kinematic distance. Column 12 
error in the kinematic distance. Column 13 effective radius (geometric mean of the semimajor 
and semiminor axis). Column 14 free- free flux. Column 15 dynamical time defined as the 
average bubble radius divided by the corrected expansion velocity. Column 16 the result of 
KDA resolution, and Column 17 the references used to resolve KDA. 
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Fig. 3. — Logarithm of fractional error in kinematic distance ^D/D versus kinematic dis- 
tance D. Open stars correspond to error solely due to the streaming motion of ionized gas, 
open circles solely due to the expected velocity spread of 15 km s^^, and filled squares the 
total error including the propagation of uncertainties in galactic longitude of the SFC 0.5 
deg, the uncertainties in the rotation curve, and the error due to spiral structure. The ob- 
jects at ~ 8 kpc with /S.D/D ~ 100 have forbidden values of v^., so the errors in Vj. have 
been taken to be At>j. = ±20 km s~^. This plot excludes all the SFCs that were manually 
assigned either to 4 kpc or 12 kpc. 
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Fig. 4. — LSR velocity plotted against galactic longitude of the centroid of the SFCs. Open 
circles denote the SFCs with kinematic distance ambiguity resolution while the filled stars 
denote those without resolution. All the SFCs with LSR velocity measurement are plotted. 
The slanted bar-like feature in |Z| < 30 and a transv erse velocity width of ±20 km s~^ corre- 



sponds to the molecular ring (see lDame et al.ll200ll . their figure 3). Star forming complexes 
in the ring account for ~ 50% of our sources with vlsr measurements, 31% of the Galactic 
free-free flux, and 35% of the free-free luminosity of the Galaxy. This is likely an underesti- 
mate since, of the 127 SFCs in the ring, 59 SFCs do not have luminosity measurement due 
to the kinematic distance ambiguity. Note that the two sources with kinematic distance am- 
biguity at I > 150° have erroneous distances; they have forbidden radial velocities (negative 
velocity expected). 
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The error in the kinematic distance is determined from the propagation of err ors in the 



deter mination of Galactic longitude (A / = 0.5°), rotation curve {A a = 5km s ^, IClemens 



19851 ). 10 km s~^ streaming motion, spread in the velocity taken as measured half- velocity, 
corrected half-velocity if determined, or 15 km s~^ if the measured half-velocity was not 
determined or was less than 3 km s~^ (RMIO considers a velocity spread less than 3 km s~^ 
as multiple measurement of the same H II region), and the uncertainties in the kinematic 
distance for SFCs associated with spiral arms or the molecular ring. In case of manual 
distance designation (for instance, some of the known SFCs near the Galactic centre; these 
SFCs are flagged in Table [1]) or for SFCs with no valid error calculation due to inaccurate 
vlsr, the error in the distance is defined as 30% of the kinematic distance, following RMIO. 
The median error contribution of each component is 1% for galactic longitude, 2% for rotation 
curve, 14% for streaming motion, 16% for the expansion velocity, and 21% for the spiral arm 
(only for the sources within Ikpc of the spiral arm). While spiral arm seems to be the most 
significant error contributor, only 72 of the 153 SFCs (all the SFCs with distances calculated 
by the rotation curve, not by manual assignment) were within 1 kpc of the spiral arm. We 
conclude that in the most general case, the streaming motion and the expansion velocity are 
the most significant error contributor. The overall median distance error is 35%. 

The RRL velocity measurements of the individual H II regions are presented in Table 
|2] where Column 1 gives the corresponding SFC number. Column 2 the name of the H II 
region. Columns 3 and 4 the galactic coordinates. Column 5 the LSR velocity, and Column 
6 the reference. 



2.3. Comparison to RMIO 

All the SFCs, including those in the 13 WMAP sources that are discussed in RMIO, 
are independently selected to ensure homogeneous classfication. As shown in Figure [H by 
resizing the 13 WMAP sources studied by RMIO, a total of 74 SFCs are identified compared 
to the previous 40: 19 original SFCs are kept, 18 are resized, reoriented, split or merged, and 
33 are completely new. For a fair comparison between RMIO and this study, matches between 
the original and the new sample are divided into two groups: those that have centroids 
overlapping other SFCs (hereafter. Type 1), and those that only have an intersection between 
ellipses (hereafter. Type 2). Schematics showing these two different types are shown in 
Figure [51 As shown in Table [3l Type 1 matched objects have, on average, distance estimates 
that differ by 16%. This corresponds to an error that is three times smaller than the typical 
error in kinematic distance in the SFC catalogue, 3.0 kpc. Using the average distance to 
our SFCs, the difference in (/,6) centroids between the two catalogs yields a difference of 
about 9pc in the bubble centroid, which is substantially smaller than the typical radius 
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of massive GMCqj. Comparing the semi-major and semi-minor axes between our sample 
and RMIO shows that the SFCs as identified in this paper are somewhat larger. The two 
samples identify either the same or physically associated objects, where some of RMIO SFCs 
are inside the SFCs found in this work. However, we note that there is a non-negligible 
fraction of un-matched complexes since we find 33 new SFCs out of 74. Type 2 matches, 
in general, are significantly more discrepant than Type 1 matches, but it is likely that the 
Type 2 matches are physically different objects. 



Results 



3.1. SFC Free-Free Luminosity 



Free-free emission is produced by photoionized gas both inside and outside the SFC shell 
walls. Assuming the origin of the ionization is associated with stars at the distance of the 
SFC, the total free-free luminosity of all the SFCs has a direct relationship to the galactic 
star formation rate (SFR). The energy p er unit volume per un it time emitted by free- free 
emission in an ionized region is given by ( iKarzas fc Latter! Il96lh : 



1/2 



(1) 



^3kb me 

where Z is the charge per ion, kb is Boltzmann's constant, c the speed of light, me the 
electron mass, T the electron temperature, and rij electron and ion densities, respectively, 
and gjf is the Gaunt factor. Following the procedure of M RIO, we adopt = rij, Z = 1, 
Te = 7000 K, typical values of H II regions, and gff = 3.28 (ISutherlandlll998l ) using 94GHz 
corresponding to W band WMAP free- free map. This corresponds to e^'^ = eonl where 
eo = 2.67 X lO^^^g cm^ s'^ Hz'^ 



The total number of ionizing photons required to keep an isotropic H II region ionized 



IS 



Qtot = j nla{H+)dV, 



(2) 



where a{H~^) = 3.57 x 10 ^'^cm ^ is the H recombination coefficient (lOsterbrock 



19891 ) and V is the volume of the ionized region. The total ionizing luminosity of a given 
H II region is then 



'"^ Since most stars are b elieved to be born in GMCs, it is likely that the SFCs presented here are embedded 
in GMCs. iMurravl (|2011l ). for instance, matches the SFCs to known GMCs to estimate the star formations 
efficiencies and lifetimes of the GMCs. 
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Q,„, = ^^^L, 1.34 X 10^6 ( photons 8"^ (3) 
eo \eTgs V 

where Ljy is the specific free-free luminosity (ergs per second per Hz), at 94GHz, of the 
given H II region. 



3.1.1. Galactic Distribution of SFC 



The galactic distribution of all the SFCs in the catalogue (left panel) and only those with 
unique distance and free-free flux measurement (right panel) are shown in Figure El We note 
that the SFC distribution in Figure |6] shows non-axisymmetric features. This is partly due 
to the paucity of RRL observations in the fourth quadrant compared to the first quadrant. 
It is also possible that the SFCs are tra cing the non- a xisyrn metric matter distribution in the 
Milky Way galaxy, previously found by iDame et al.l (119871 ) . 



Past studies have shown that most star formation in the Galaxy occurs in molecular 
clouds, and the Galactic distribution of molecular hy drogen falls off rapidly beyond the solar 
circle (IDame et al.ll2001t iBinney fc Merrifieldl 119981 ). Thus we should expect that the star 
formation falls off with distance from the galactic center, at least beyond the molecular ring. 
This behavior is apparent in Figure H] where ~ 50% of all the SFCs with LSR velocity mea- 
surement, and 35% of the free-free lumi nosity, are associated with the molecular ring, which 
is at Galacto-centric radius ~ 3.6 kpc ( IBinney &: Merrifieldl Il998l ). As expected. Figure [6] 
shows that there are significantly fewer SFCs outside the solar circle (hereafter, OSC) than 
those inside the solar circle (hereafter, SC). 

We have compiled the luminosity function of star forming complexes inside and outside 
the solar circle in Figure [71 A Kolmogorov-Smirnov (KS) test shows that the normalized 
luminosity functions do not differ at ~ 100%. 

While the shapes of the luminosity functions inside and outside the solar circle are sim- 
ilar, the total free-free luminosity of sources inside the solar circle is 1.0 x 10^^ erg s~^ Hz~^; 
that of sources outside the solar circle is only 0.2 x 10^^ erg s~^ Hz~^. Furthermore, the 
distances to sources outside the solar circle may well be overestimated by a factor of 1.5, as 
mentioned in § 12. 2^ so that the corrected OSC cumulative luminosity might be as little as ~ 
9% of the cumulative luminosity of sources inside the solar circle. 



3.1.2. A Possible Distance Bias 



Figure [61 shows that the surface density of star forming complexes is higher near the 
sun than at large heliocentric distances, suggesting that we may be missing sources at large 
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Fig. 6. — Left: Galactic distribution of all the SFCs with kinematic distances. The dashed 
circle is the solar circle and its intersection with the two dashed lines is the position of the 
sun. Near and far distances of a single source are connected with dotted lines. A solid circle 
denotes SFCs with a unique distance measurement or KDA resolution. Right: Galactic 
distribution of all the SFCs with free-free luminosity determination. Colours represent the 
free-free luminosity in erg s~^ Hz~^ on a logarithmic scale. 
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distances and so underestimating the Galactic free-free luminosity. To see if this is the case, 
we examine the luminosity function and the cumulative luminosity of sources on the near 
and far side of the Galaxy. By the near side of the Galaxy we mean sources with y > in 
Figure [6l sources with y < are then on the far side of the Galaxy. The luminosity function 
and cumulative luminosity as a function of Q (for near-side star forming complexes, far-side 
complexes, all complexes inside the solar circle, and all the SFCs) are compared in Figure [7] 
and Figure El respectively. 

Figure [7] shows that the total number of far-side sources is 38, less than half the total 
number of those on the near-side (84), consistent with the impression given by figure |6l 
Furthermore, Figure [7] shows a turnover at Q ~ 10^*^ s^^ — the break is seen more clearly in 
the differential luminosity function, to be discussed below. In fact, we find no sources with 
Q < 10^^ s on the far side of the Galaxy. This demonstrates that we have a selection bias 
against l ow luminosity (Q < 10^s~^) sources at large distances, a result similar to that 



found by [Murray fc RahmanI (120101 ) 



Despite this. Figure [H] shows that the total ionizing luminosity of sources on the far-side 
is ~ 31% greater than that of near-side sources. In fact. Figure [8] shows that the 24 most 
luminous sources (with dust-uncorrected Q > 1.3 x 10^^ s ^^), out of 280 total sources (~ 9%) 
account for a cumulative (dust-uncorrected) Q ^ 1.5 x 10^^ s~^, or 75% of the luminosity 
from sources with known distances. These sources account for half of the total free-free 
luminosity of the Galaxy, as determined below. So while it is true that our selection has 
a bias against distant, low luminosity sources, low luminosity sources make an insignificant 
contribution to the total free-free emission of the galaxy. 

Figure [3 shows the same data as in figure [HI but this time integrating from the minimum 
Q to large Q, making it easier to measure the slope of the cumulative luminosity function 

The data can be described by two power laws, one for Q > 6.7 x 10^° s~^, and another for 
Q < 4.2 X 10^°s~^ (there is a hint o f a thi rd break at Q < 10^ " s~^ b ut it is not very clear). 



Following iMcKee fc Williamsl (119971 ) and [Murray fc Rahman! (120101 ). we define the power 



law exponent F by 

oc Q-^, (5) 

and Q oc , where F = 1 — F. The red line in figure [9] is a least squares fit to the data above 
Q ^ 6.7 X 10^^ s -\ with f ^ 0.7. The green line shows that f ^ 1.5 for Q < 4.2 x 10^° s -\ 
In other words, for Q < 4.2 x 1050 s~^, F = —0.5, while for larger Q we find F = 0.3. 

As a check, figure [TO] shows the differential luminosity function dN/dQ. Again, we find 
two different slopes, defined as (A logrfA^)/(A logQ) = — F. for Q < 4.2 x 10^^ s~^, we find 
F !^ —0.5 and F 0.7 for larger values of Q. These exponents are in reasonable agreement 
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Fig. 7. — The cumulative ionizing luminosity function (dust-uncorrected) of SFCs with 
unique distance and flux measurement. Open circles denote SFCs inside the solar circle on 
our side of the Galaxy > in figure [6l labeled NEARSFC): solar circle SFCs on the far 
side (y ^ 0) are denoted by open stars (FARSFC). Filled squares denote all objects inside 
the solar circle (SC). The open diamonds (ALL) denote all the SFCs, both inside and outside 
the solar circle, with the distance correction factor of 1.5 as discussed in § l2.2.1l applied when 
calculating Q. 
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Fig. 8. — Cumulative ionizing luminosity versus ionizing luminosity (dust-uncorrected). The 
legend is the same as Figure [71 Note that the cumulative luminosity of sources inside the 
solar circle, but on the far side of the galaxy from the sun, exceeds that of sources in side 
the solar circle but on the near side of the galaxy by ~ 31%. The 24 most luminous objects 
(those with Q > 1.3 x 10^^ s"^) produce a total Q = 1.5 x 10^^ s~^, half the (non-dust 
corrected) ionizing flux of the entire Galaxy. 
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with the corresponding values obtained from figure M While F of smaller Q's are in perfect 
agreement, the exponents for larger values of Q agree within ~ ±0.4 precision. We note that 
half of the bins in the high Q regime in [10] are on the order of 1 or 2 and suffer from poor 
sampling. In fact, we find the exponents to change by a factor of two within two decades of 
binsizes. We argue that the two luminosity functions are in agreement within the precision 
of the sampling. 

The slope F = 0.3 — 0.7 of the high luminosity tail is similar to the values, ranging frorn 
0.41 to 1.3, with an average F = 0.75, found in other galaxies by iKennicutt et al.l (119891 ) 



and McKee fc WiUiams (1997). Whi 



e the slope we find is shallower than the mean value 



found bv I Gonzalez Delgado &: PerezI (1997), i t is in agreement with their slope for later type 



galaxies (see iGonzalez Delgado Sz Perezlll997l . their Figure 3). The value of F we find for the 
Milky Way is consistent with our result that 24 sources are responsible for half the ionizing 
luminosity of the Galaxy. 

We conclude that our sample is biased against finding low luminosity sources at large 
heliocentric distances (larger than a few kpc), but that this bias has only a minor effect on 
our estimate of the Galactic free-free or ionizing luminosity. 



3.1.3. Accounting for sources with unknown distances 

The 168 SFCs that have both unique distance and free-free flux measurement only 
account for 58% of the total Galactic free-free flux. The 280 SFCs that have free-free flux 
measurements (but not necessarily unique distance measurements) account for 74% of the 
total Galactic free-free flux (see §2.2.21 above). To estimate the luminosities of the 112 SFCs 
with known flux but unknown distances, we perform a Monte Carlo (MC) simulation. To 
do so, we need a model for the galactocentric radius distribution. We assume that the star 
formation is distributed in a manner similar to that fo und for molecular gas, as shown, for 



example, in Figure 9.19 of iBinney fc Merrifleldl ( 1l998l ). In particular, we model the SFC 



distribution as a Gaussian with mean radius /i = 5 kpc and dispersion a = 5. 7 kpc. These 
parameters, and the dispersion in particular, are found by requiring that the number of 
SFCs inside and outside the solar circle match the observed numbers for our sample of 168 
sources. For each SFC without distance, a galactocentric radius is randomly sampled from 
the gaussian distribution reflecting the variations of the molecular gas. If the MC-generated 
radius is not allowed at the measured Galactic longitude of the SFC, we repeat the process 
until we find an allowed radius. 

Fixing the Galactic radius and longitude of the MC source leaves us with a Heliocentric 
distance ambiguity if the source is inside the solar circle. To resolve this ambiguity, we apply 
three criteria. First, the free- free luminosity of the SFC must be greater than the minimum 
luminosity detectable at that distance. We estimate this minimum luminosity from the data; 
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Fig. 9. — The cumulative ionizing luminosity Q{< Q) versus ionizing luminosity Q. The 24 
most luminous objects (those with Q > 1.3 x 10^^ s~^) produce a total Q = 1.5 x 10^^ s~^, 
half the (non-dust corrected) ionizing flux of the entire Galaxy. 
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Figure [TT] shows the minimum detected luminosity as a function of Hehocentric distance for 
the 166 sources with known distances. We fit a power law to the data, and use this fit to 
determine the minimum detectable luminosity as a function of distance. We find that this 
minimum detectable luminosity scales roughly as (the actual value of the exponent is 
1.96) suggesting that there is in fact a roughly fixed flux limit. We refer to the requirement 
that our MC sources exceed this minimum luminosity as the sensitivity check. 

The second criterion we use to resolve the distance ambiguity relies on our knowledge of 
Galactic structure; we note that SFCs have a higher chance of being associated with spiral 
arms than of being between arms. If the distance ambiguity is not resolved by the sensitivity 
check, then the d i stanc e that is within 1 kpc of an arm in the galactic spiral arm model of 



Taylor fc Corded ( 119931 ) is chosen; we call this the spiral arm check. 



Finally, we expect most, if not all, SFCs to lie within the Galactic plane. We assert 
that the height of the SFC above the plane, z, calculated using near distance {znear) must be 
smaller than 60pc (scale height of the Galactic molecular disk), and the z calculated using 
far distance (-Z/ar) be smaller than 200pc. If the ambiguity remains after sensitivity and 
spiral arm check, and Zmar < 60pc and Zfar > 200pc, we adopt the near distance; we call 
this height check. These parameters are found by requiring that the number of SFCs on 
our side and the other side of the Galaxy inside the solar circle match the observed number 
for our sample of 168 sources. While we expect most SFCs to lie within the scale height 
of the molecular disk, 60pc, above the plane, especially for those SFCs associated with the 
molecular ring, we leave the possibility of outliers. In fact, 16 SFCs out of the 168 SFCs 
with unique distance measurement have z > 60 pc. Of these 16 SFCs, 6 of them are likely 
associated with the molecular ring, and only one of these 6 SFCs is at z > 100 pc (but its 
Galactic latitude is 1.8°). We also note that the height has a mean uncertainty of ~ 35% 
due to the uncertainties in the distance derived from RRL velocities. 

If both distances pass or fail the height check, we choose the distance assuming a uniform 
probability distribution between the two distances, as we have no other reason to favour one 
side over the other. 

To check our model, we perform MC simulations with 168 SFCs having galactocen- 
tric radii drawn from the normal distribution described above. We then assign a galactic 
longitude chosen randomly from the allowed range governed by the galactocentric radius. 
Galactic latitudes are drawn from the normal distribution with fi = 0.11° and a = 0.90° 
where the parameters are the results of the gaussian fit of the full sample of 280 SFCs. 
Finally, we assign a free-free fiux randomly from the catalogue of fiuxes from the 168 SFCs 
with distances. Near and far distances of the SFCs within the solar circle are then assigned 
based on the three checks described above. 



Repeating this process 5000 times result in, on average, 91 SFCs are on the near side 
of the Galaxy, and 30 on the far side, all inside the solar circle. This leaves 47 SFCs outside 
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Fig. 10. — The differential luminosity function of SFCs . The y-axis gives the number of SFCs 
in logQ = 0.176 bins. Filled (red) circles are the observed SFC sample (166 objects) while 
open squares represent the MC simulated SFCs. A KS test shows that the two samples 
do not disagree at the 20% confidence level. The open circles with error bars represent 
the differential luminosity function generated by combining the 168 objects with known 
luminosities with a MC model for the 112 objects lacking measured distances. 
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Fig. 11. — The minimum luminosity of observed SFCs plotted against the (binned) kinematic 
distance. The slope of the fitted hne is a 1.96. 
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the solar circle. This is similar to the observed numbers, by construction. A KS test shows 
that the simulated and observed differential luminosity functions do not disagree at the 20% 
confidence level. Figure HO] shows the differential luminosity function of the 168 observed 
SFCs with known luminosities (filled circles), as well as that generated by our MC model 
(open squares). 

Repeating the MC simulation 5000 times for the remaining 112 SFCs which have free- 
free flux measurements but lack distance determinations, the model predicts 94 SFCs inside 
the solar circle, with 80 on our side of the Galaxy and 14 on the far side. The free-free 
luminosity of the Galaxy ("All"), of sources inside the solar circle ("SC") and on the near 
and far side (but inside the solar circle) predicted by the MC simulations, are all summarized 
in Table HI which also lists the observed luminosities. The table heading "Blind-MC" refers 
to the MC sample of 168 SFCs with randomly assigned Galactocentric radii (chosen from the 
normal distribution described above), longitudes chosen randomly from the allowed range, 
latitudes chosen from the normal distribution described above, and free-free fluxes randomly 
assigned from the measured fluxes of the 166 sources. The table heading "Data-MC" refers 
to the combination of the 166 SFCs with measured luminosities, supplemented by the MC 
simulation of 112 SFCs using the observed Galactic longitudes and free- free fluxes. 

The MC estimate of the free-free luminosity (Blind-MC) agrees with the observed free- 
free luminosity (both for 168 SFCs) at the one sigma level, as shown in Table |H We conclude 
that our MC model is reliable at the la level, and henceforth use the luminosities from the 
third column of table H] (for all 280 SFCs) as our best estimate for the free-free luminosities 
from the corresponding regions of the Galaxy. 

3.2. Luminosity and SFR 

The Galactic free-free luminosity due to star forming complexes is given by the sum of 
the free-free luminosities of the 168 sources with measured fluxes and distances, supplemented 
by the Monte Carlo estimates for the luminosities of the 114 SFCs with measured fluxes but 
not distances. We find 

L^ ..pc = 1.3 ± 0.2 X 10^^ erg s"^ Hz"^ (6) 

at W band, where the error estimate includes the error in the kinematic distance as described 
in § l2.2.5l added in quadrature to the la error of the Monte Carlo simulation. The luminosity 
of objects outside the solar circle are reduced by the factor of 1.5 described in § 12.2.11 without 
this correction, L^^sfc is larger than that given in eqn. |6]by ~ 10%. 

We note that our estimate of the ionizing luminosity of the SFC associated with the 
Arches cluster (G359 SFCO) is 1.8 x 10^^ s~^ which is about 4 times less than expected 
assuming Arches cluster has 150 03 main sequence stars (6.5 x 10^^ s^^). This is likely due 
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to an underestimate of the free-free flux assigned to SFCO. This is a warning that, as the 
Galactic Centre is a highly confused region, the SFC parameters in that region should be 
taken with caution. 

The total free-free flux in 280 SFCs is F^,sfc = 34234.8 Jy, or 74% of the galactic 
free-free flux Fj^^gai = 46177.6 Jy. This result is consistent with Figured! which shows that 
most of the free-free flux is concentrated along the Galactic plane where we find SFCs. We 
believe that the remaining 26% of the flux arises from the diffuse component of the free-free 
emission leaking out of the WMAP sources. If so, then the unaccounted for flux is emitted, 
on average, at very nearly the same distance as that emitted by the SFCs. Hence we estimate 
the total Galactic free-free luminosity (due to both SFCs and to diffuse emission) to be 

Lu = L,,SFC + ^irD^F.^gai - F,,sfc) = 1-6 ± 0.3 x 10^^ erg s'^ Hz-\ (7) 

where 

Df = ^^f' = 4.7± 1.8kpc (8) 

is the flux weighted mean distance to the 168 SFCs for which we have both unique distance 
and flux estimates, supplemented by the Monte Carlo estimates for the luminosities of the 
114 SFCs with measured fluxes but not distances (the value changes to 4.8 kpc without the 
Monte Carlo estimates). If we assume the distance is 8 kpc, as expected if most SFCs are 
in the 3 kpc molecular ring), Li, changes to 2.2 ± 0.3 x 10^^ erg s~^ Hz~^, 38% higher than 
the quoted value. Without the correction for the luminosity of the objects outside the solar 
circle, increases by ~ 8%. 

The error estimate in equation [7] comes from propagating the errors in the distance esti- 
mates. The error in the mean distance contributes an error of 25%, the error in the distance 
of individual SFCs contributes another 8%, while the la error in the MC simulation accounts 
for 1%. In § 12.2.51 we showed that the most significant error in the kinematic distances were 
due to the expansion velocity of SFCs. We conclude that the most significant source of 
uncertainty in our estimate of the Galactic free-free luminosity is the expansion velocity of 
SFCs, combined with the uncertainty in the mean distance of the sources responsible for the 
diffuse flux. 

Using EquationEl the total ionizing luminosit y, before the correction for dust absorption, 
is Q = 2.1 ± 0.4 X 10^^ photons s~^. Following iMcKee &: Williams! ( 119971 ). we correct for 



the effect of ionizing photons' absorption by dust grains by multiplying by 1.37: 

Q = 3.0 ± 0.5 X 10^2 photons s"^ (9) 

compared to the estimate of MRIO, Q = 3.2 ± 0.5 x 10^'^ s To calculate the star formation 
rate of the Milky Way, we use the conversion from ionizing flux to stellar mass given by 
MRIO: 

M,/Q = 4.1 X lO-^^Moyr-^s. (10) 
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The star formation rate of the Milky Way is then 

= 1.2 ± 0.2 Mq yr~\ 



(11) 



which differs from the estimate of MRIO by ~ 8%. As MRIO note, however, the estimate 
given above is highly dependent on the choice of the IMF, in par ticular the high mass end 
slope (denoted by F). We follow MRIO and assume F = 1.35, the ISalpeterl (Il955l ) value. 



This estimate for the sta r formation rate also agrees, within one standard deviation 



with the SFR determined by iRobitaille fc WhitnevI feoiol ). 0.68 - 1.45 M© yr"^ These 



< 



authors determine the star formation rate by counting intermediate to high-mass (SMq 
M < 25M0) pre- main sequence stars, in contrast to our use of a tracer of > 40Mq stars, 
which are responsible for the bulk of the ionizing flux produced by young star clusters. 

We note that SFCs inside the solar circle account for 80% of the total galactic free-free 
luminosity, i.e., they are responsible for roughly 80% of the star formation activity in the 
Milky Way galaxy. 



3.3. The relationship between free-free and 8/im emission 

With a more complete sample, we re-examine the relationship between WMAP free-free 
flux and 8yum from RMIO as shown in the top panel of Figure [121 

Fs,m oc (12) 

where Fg^m is the aperture-integrated 8yU.m flux and Fy is the aperture-integrated 94 GHz 
free-free emission from the WMAP free-free foreground emission map. The MSX 8/im flux 
is calibrated to GLIMPSE by an empirically determined relationship: 

Fgl(J?/) = 5.19436 ± 0.09625 FA/5x(^y) - 5921.22 ± 2587.89. (13) 

The calibration is performed by first retrieving the total fiux of overlapping regions between 
GLIMPSE and MSX images then fitting a linear line between the two fiuxes. The correlation 
coefficient of the fit is 1.0. 



Calzetti et al.l (120071 ) find the power law in Equation IT^ to be 0.94 ± 0.02 in the studies 
of extragalactic star formation, with the caveat that the analysis should be limited to the 
regions of star formation with a narrow spread in metallicity around the solar value. 



As ICalzetti et al.l (120071 ) note, 8/im emission traces the stellar light absorbed and re- 
emitted mostly by polycyclic aromatic hydrocarbon (PAH). Since both free-free and 8/im 
emission are tracers of star formation, a linear relationship is expectecj^. 



^Due to a coding error, RMIO reported that the 8/xm flux scaled as the square of the free-free flux. 
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The flux relationship for SFCs follows a similar power law, as shown in the bottom panel 
of Figure [121 

Fs,n. oc Fr^'-'^'. (14) 

We have also found a tight correlation between the 8/im emission and the radio recom- 
bination line flux in our SFCs: the correlation coefficient is 0.7. 

Finally, we note that the SFC 8fim flux exhibits a strong correlation with the projected 
size of the SFCs. Nearby sources of a fixed luminosity and physical size will appear both 
larger and have larger fiuxes than equivalent sources at larger distances. If this were respon- 
sible for the correlation, one would expect that the fiux would scale with the area on the 
sky (the surface brightness would be constant), whereas we see that the fiux scales with the 
square root of the area. 

A possible physical reason for such a c orrelation is the destruction of PAH molecules 



near the ionizing source (IWatson et al.ll2008l ). reducing the 8/im emission in compact SFCs. 
MRIO also finds that the 8/im fiux is dominated by emission at large projected radii rather 
than by the small, bright areas toward individual H II regions. 

The sum of the 8/im emission inside SFCs, i.e., inside the bubble walls, only accounts 
for 33% of the 8/tm fiux within the (more extended) WMAP sources. This suggests that the 
far UV (FUV) photons that excite PAH emission travel well past the bubble walls defining 
the SFCs. However, there may be sources other than central clusters that are responsible for 
the PAH emission in the diffuse region outside of the SFCs. One indication of this is that 
we find larger scatter in the 8/im fiux when we sum over a larger aperture in the GLIMPSE 
images. It is possible that nearby old B stars (with no association to SFCs) are emitting 
FUV photons. 



3.4. SFC Expansion and Molecular Gas Turbulence 

RMIO suggest that the primary driver of turbulence in Milky Way molecular clouds is 
the expanding bubbles seen inside SFCs. We test this idea using the more complete sample 
presented here. 

Equation 3 of RMIO is used to determine the mechanical luminosity in the expansion 
of each SFC: 

Lmech = 2^oAi7fr, (15) 



where Eq = 1OO±88M0 pc"^ (IHever et al.ll2009r i is the surface density of a CMC (while their 
result shows an LTE-derived median surface density as 42 ± 37 Mq pc~^, they consider 
this a lower limit and that the real value can change by a factor of 2-3 due to their neglect of 
abundance variations in the outer envelope of the clouds), Avc corrected half-spread velocity. 
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Fig. 12.— Top: Log-log plot of 8//m vs. free-free flux of the WMAP sources with GLIMPSE 

or MSX coverage. The free-free flux of the WMAP sources is defined as the maximum 
between the isophotal and aperture-integrated 94 GHz free-free emission. The slope is found 
to be 0.73 ± 0.04 with a correlation coefficient 0.93 and reduced = 4.87 which shows that 
the two emissions arc strongly correlated by approximately a linear relation. The error in 
the -Fg/im is defined as the difference between the ffux within SExtractor-determincd aperture 
and the flux within the resized and reoriented aperture. Bottom: Log-log plot of 8//.m vs. 
free-free flux of the SFGs within the WMAP sources with GLIMPSE or MSX coverage. The 
slope is found to be 0.87 ± 0.07 with a correlation coefficient 0.71 and reduced — 3.44 
which shows a relatively weak correlation between the two emissions. The error is determined 
by taking the standard deviation of the total flux within apertures ranging from 0.8 to 2 
times the original SFC radius; the mean error among the SFGs, 88%, is used. 
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and r the geometric mean radius. For the SFCs without a measurement of the corrected 
expansion velocity, we use a mean expansion velocity of 11.7 ± 5.9 km s~^. 

Limiting the analysis to SFCs inside the solar circle and correcting for the SFCs that 
do not have a KDA resolution as in § I3.1.2[ the total mechanical luminosity being injected 
into the interstellar medium (ISM) inside the solar circle is 



L^ech = 1.0 ± 0.2 X 10 



39 



( 



lOOMopc- 



erg s 



(16) 



where the error estimate results from propagation of the errors in the geometric radius, 
the expansion velocity, and the surface density. The contributions are 5% (radius), 17% 
(expansion velocity), 12% (surface density), and 1% from the MC simulation in § 13.1.21 
Once again, the most significant contributor is the error in the expansion velocity. 

This bubble kinetic luminosity should be compared to that required to maintain the 
velocity dispersion in the molecular gas within the solar circle: 



■'turb „ 



lMf2 



(17) 



where M is the total molecular gas mass inside the solar circl e , v rn olecular gas velocity, and 
to the turbulence dissipation time scale defined by lMac Lowl (119991 ): 



1.0- 



(1^ 



where C is the driving wavelength. We take the molecular gas velocity as the rms velocity 
and choose 2h, twice the scale height of the molecular disk, as the turbulence driving scale 
since it is the largest mode wavelength of a marginally gravitationally stable differentially 
rotating disk (iGoldreich fc Lynden-Belllll965l ). 



Using this, the turbulent luminosity is 



IMw^ 



-'turb 



2 2h 



(19) 



We employ h = 60 dc (iMalhotral I1994J 1. M = 1.0 x 10^ Mq, and v = V2 In 2amoi with 
Cmoi = 7 ± 1 km s~^ ( Malhotral I1994I ) . The estimated solar circle molecular gas turbulent 
luminosity is then 

3 



■'turb 



1.5 ± 0.6 X 10 



39 



7 km s 



60pc 



erg s 



(20) 



The kinetic luminosity due to the expansion of the bubbles inside the SFCs can ac- 
count for 65% of the solar circle molecular gas turbulent luminosity. This suggests that the 
expansion of the bubbles inside SFCs is a major driver of turbulent motions of molecular 
gas. 
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4. DISCUSSION 

We note that the bubble expansion velocity is larger than the Galaxy-wide average 
turbulent velocity, while the average bubble radius is smaller than the scale height of the 
molecular disk. It follows that the bubble time scale is significantly shorter than the turbulent 
decay time. Thus the total energy contained in the turbulent motions is significantly larger 
than that contained in the bubbles. It is also clear that the mass swept up in the bubble 
walls is significantly smaller than the total molecular gas mass inside the solar circle. 

Estimates of the virial par ameters of the mo st massive GMCs suggests that they are 



gravitationally bound objects (IHeyer et al.l 120091 ) . yet the turbulence that supports them 



decays on a crossing time. Since the clouds are known not to collapse (such a collapse 
would lead to very high star formation rates) then either some way to sustain the turbulence 
is needed, or the GMCs must be disrupted. The fact that we find expanding bubbles in 
essentially all the rapidly star forming GMCs in the Milky Way strongly suggests that the 
clouds are disrupted rather than maintained in a quasi-steady state by injection of turbulent 
energy. 

The bubble expansion can be driven by a number of possible mechanisms, including 
radiation pressure from star clusters in the bubble, warm i onized gas pressure , and hot (x- 



ray emitting) gas pressure (IHarper-Clark fc Murrayl l2009l : iMurray et al.l l2010l : iLopez et al. 



20 111 : iPellegrini et al.ll201ll ). It is possible that some of the low-luminosity SFCs may have 
low ionizing luminosities because they are > 4 Myr old, and hence may contain supernova 
remnants. In fact, a few SFCs are observed to contain known remnants and show bright 
emission in the WMAP synchrotron map. A uick visual inspection of WMAP W band 
synchrotron emission map shows that only about 23 SFCs (8%) overlap bright synchrotron 
emission map. This suggests that supernovae do not disrupt GMCs. It also suggests, less 
directly, that supernovae may not be the dominant non-gravitational source of stirring in 
the ISM of the Milky Way. A more rigorous search for synchrotron emission in the faint, 
possibly older SFCs and WMAP sources would be very instructive. 

When we say that bubble expansion is a major driver of turbulent motions, we have 
in mind the following picture. Imagine starting with a turbulently supported gas disk with 
scale height /i, turbulent velocity cr, and a gas surface density S. Imagine further that the 
Toomre Q = Vcc/lnGUr) parameter is near unity, so that the disk is marginally unstable 
due to self-gravity (fc is the circular velocity of the galaxy, and r is the galactocentric 
radius). Over time the turbulence decays, so that a decreases, as does Q. This triggers 
large scale instabilities, which drive turbulence on the scale h. However, some fraction of 
the gas collapses to form gravitationally bound objects, i.e., giant molecular clouds. These 
clouds contain subregions which are also gravitationally bound (often called clumps in the 
literature). These clumps are where the massive star clusters we study are formed. The 
star clusters emit radiation and stellar winds, which both act to inflate the bubbles we see. 
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The bubbles expand rapidly, eventually disrupting the host GMC, and lifting gas above and 
below the galactic disk. Finally, the material falls back towards the disk plane, once again 
providing a source of large scale turbulence via gravitational instability. 



5. Conclusion 

We have analyzed 72 of the 88 Galactic WMAP sources and found 280 star forming 
complexes. Using 8/im fluxes, we assign free-free fluxes to all 280 SFCs. Of these 280 sources, 
168 have a unique kinematic distance measurement. The 112 SFCs which have free- free flux 
estimates but which lack a unique kinematic distance measurement are accounted for by 
assuming a model of Galactic spiral arm, molecular hydrogen distribution and a detection 
sensitivity as a function of galactocentric radius. These 280 SFCs account for 74% of the 
free-free flux of the Galaxy. The remaining 26% of the free-free flux is shared between the 
flux associated with the sixteen WMAP sources lacking complete 8/im coverage or bubbles, 
and the flux from outside the WMAP sources. 

While our sample is biased against low luminosity sources at large heliocentric distances, 
we find that low luminosity sources have only a minor contribution to the total Galactic free- 
free luminosity. In fact we show the most luminous 9% of all SFCs account for half the total 
Galactic free-free luminosity. The SFC luminosity function is described by two power-laws, 
with the slope in the high-luminosity tail consistent with values found in the literature on 
extragalactic H II luminosity functions. 

With the SFC sample in hand, we extend the results of MRIO and RMIO. The total 
galactic free-free luminosity is measured to be 1.6 ± 0.3 x lO^^erg s~^ Hz~^, and the corre- 
sponding ionizing luminosity is 2.9±0.5 x 10^^ photons s We find a Galactic st a r form ation 



rate of 1.2 ± O.2M0 yr -\ in agreement with MRIO and iRobitaille fc WhitnevI (120 lOf ). We 
find a nearly linear relationship between 8/im and free- free flux of the WMAP sources and 
the SFCs, consistent with extragalactic studies. Solar circle SFCs are found to be responsible 
for 80% of the total galactic free-free luminosity, and their kinetic energy injected per unit 
time into the ISM via expansion of the bubbles is equal to 65% of the luminosity required to 
maintain the velocity dispersion seen in the solar circle molecular gas. This result supports 
the suggestion of RMIO that the expansion of the bubbles due to central massive clusters 
inside the SFCs is a major driver of the turbulent motion of the molecular gas in the inner 
Milky Way. 
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Table 1. Star Forming Complex Parameters 



No. 


WMAP 
NAME 


I 

(dcg) 


b 

(dcg) 


smaj 
(arcmin) 


smin 
(arcmin) 


PA 
(deg) 


vlsr 
(km s"-*-) 


A Dm 

(km s"-"-) 


(km s"-*-) 


D 
(kpc) 


(kpc) 


<R> 
(pc) 


U 

(Jy) 


^dyn 

(Myr) 


KDA 
N/F 


Ref. 





G359 


0.129 


0.044 


4.6 


3.4 


-30 


-21.7'^ 




25.5<= 


8.0 


2.4 


9.2 


127 


0.4 






1 


G359 


0.165 


-0.060 


2.3 


1.7 


15 


34.0^* 




48. 7^* 


8.0 


2.4 


4.6 


34 


0.1 






2 


G359 


0.285 


-0.501 


20.5 


11.8 


15 


20.0 


3.2 




9.0 8.0 


1.5 1.5 


40 fi 36 


466 


3.5 3.1 






3 


G359 


0.502 


-0.056 


2.5 


1.6 


35 


47.1 






8.9, 8.1 


1.5, 1.5 


5.2, 4.8 


31 








4 


G359 


0.510 


0.183 


1.4 


1.0 

















4 








5 


G359 


0.665 


0.644 


5.2 


3.9 


60 


3.7 






11.6, 5.4 


7.5, 7.5 


15.2, 7.1 


25 








6 


G359 


0.829 


0.207 


2.4 


1.3 


-65 


9.7 






10.5, 6.5 


2.6, 2.6 


5.4, 3.4 


6 








7 


G359 


1.137 


-0.087 


3.1 


2.1 


15 


-21.7 












21 








8 


G359 


1.341 


0.118 


2.7 


2.4 


50 


-12.0 






12.0 


3.6 


8.9 


13 








9 


G3 


2.268 


0.234 


3.6 


2.4 


-10 


4.9 






14.0, 4.0 


1.2, 1.2 


12.0, 3.4 


6 








10 


G3 


2.515 


-0.032 


1.7 


1.3 


40 


8.3 






13.0, 4.0 


3.7, 3.7 


5.5, 1.7 


1 








11 


G3 


2.602 


0.162 


2.4 


1.6 


-70 


102.4 






9.2, 7.8 


0.1, 0.1 


5.3, 4.5 


3 








12 


G3 


2.866 


-0.032 


5.1 


3.2 


-40 


-1.8 






12.0 


3.6 


14.1 


14 








13 


G3 


3.297 


-0.040 


7.2 


4.8 


-40 


6.0 


4.0 




12.0, 2.6 


3.6, 3.6 


20.6, 4.5 


32 


1.8, 0.4 






14 


G3 


3.656 


-0.113 


1.5 


1.0 


-40 


4.6 






12.0, 1.9 


3.6, 3.6 


4.4, 0.7 


1 









00 



''This SFC selection is motivated by ICotera et al. I 1I2OO5I ') and iLang et al. 

^LSR is determined by taking the average between filaments El, E2 and the Arches cluster 

velocity llCotera et al.ll2005 '). 

'^The expansion velocity is determined by taking the mean FWHM velocity among E2, Wl, W2 filaments jCotera et al.ll2005l : iLang et al.ll200llV 
Velocities are taken from lLane et al" I lll997l) . 



References. — (l) [Anderson fc Banial tOO^): ( 2) iBusfield et al] 1I2OO6I) : (3) ISewilo et all (Hooi); (4) iFish et al.l tOOA) : (5) iKolpak et all tOOA) : (6) IPalagi et all ^9^; (7) 
ICaswell fc HavnesI \l98l\ ): (8) iGillespie et all lll977l ) 



Note. — The KDA Ref. indicates the literatures used to resolve kinematic distance ambiguity. In cases where no reference is given, the kinematic distance is either unique, 
undetermined, or we present both the near and far distances. 



Note. — (This table is available in its entirety in a machine-readble form in the online journal. A portion is shown here for guidance regarding its form and content.) 



Table 2. HII Region RRL Velocities 



SFC 


Name 




I 


b 


VLSR 


Ref. 


2 


[KC97c] GOOO.3-00.5 





2840 


-0.4781 


20.0 


5 


2 


[KC97c] GOOO.4-00.5 





3940 


-0.5400 


24.0 


5 


2 


[KC97c] GOOO.5-00.7 





4890 


-0.6681 


17.5 


5 


2 


MSX6C GOOO.5799-00.6302 





5718 


-0.6279 


20.0 


5 


3 


GAL 000.49-00.05 





4905 


-0.0600 


45.8 


13 


3 


[WWB83] GOOO.51-0.04 





5137 


-0.0394 


47.1 


13 


3 


[KC97c] GOOO.5-00.1 





5180 


-0.0647 


47.1 


5 


5 


[LPH96] 000.640+0.623 





6403 


0.6239 


3.7 


7 


6 


[KC97c] GOOO.8+00.2 





8290 


0.1934 


9.7 


5 


7 


[KC97c] GOOl. 1-00.1 


1 


1488 


-0.0618 


-21.7 


5 


8 


GAL 001.32-f00.09 


1 


3228 


0.0860 


-12.0 


1 


9 


[KC97c] G002.3-K00.2 


2 


3031 


0.2430 


4.9 


6 


10 


[LPH96] 002.510-0.028 


2 


5099 


-0.0281 


8.3 


7 


11 


[KC97c] G002.6-f00.1 


2 


6110 


0.1351 


102.4 


6 


11 


IRAS 17480-2636 


2 


6130 


0.1340 


102.4 


6 



Caswell fc HavnesI (Il987h : f2iHeiles et al 
()200l[ k (4jKuchar & Clark! (Il99 7|) 87G& 



Not e. — References: 
( 1996h: C3lKimfcKool , , _ ,.. _ , _ 

(5lKuchar fc ClarkI (Il997l) PM N: felLockmau (19891: f7lLockinan et al 



(1996 



(10 
(13 



Wilson et al.l 



Wink et al 




d): fsjPaladini et all (I2OO3I): f9lQuireza et all tooW : 



197C 
(1982 



q ilCaswelll (I1972I): (12jDieter_ 



(|2007h : (le iRusseil et all (12005 ) 



1967): 



(14' Peeters et all (|2002l) : fl5 1Arava et al 



(This table is available in its entirety in a machine-readble form in 
the online journal. A portion is shown here for guidance regarding its 
form and content.) 



Table 3. Star Forming Complex Comparison 





I 


b 


dsmaj 


dsmin 


dVLSR 


dAVm 


dAvc 






(deg) 


(deg) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


Type 1 


-0.04 ± 0.14 


0.10 ± 0.16 


19.77 ± 23.11 


17.51 ± 23.25 


14.39 ± 30.35 


18.96 ± 20.29 


22.50 ± 13.93 


16.09 ± 37.76 


Type 2 


-0.19 ± 0.30 


-0.13 ± 0.21 


327.33 ± 338.95 


328.95 ± 327.11 


6.58 ± 7.26 


75.69 ± 27.78 


61.01 ± 33.67 


41.77 ± 34.59 


All 


-0.07 ± 0.17 


0.04 ± 0.14 


67.79 ± 129.19 


70.89 ± 156.38 


29.14 ± 57.76 


28.32 ± 21.19 


27.33 ± 19.43 


22.13 ± 32.43 



Table 4. Monte Carlo Simulation Result 





Data 


Blind-MC 


Data-MC 


LaLL,27 


1.15 ± 0.15 


1.13 ± 0.28 


1.27 ± 0.02 


LsC,27 


0.96 ± 0.12 


0.96 ± 0.26 


1.16 ± 0.02 


LnEAR,27 


0.41 ± 0.09 


0.40 ± 0.12 


0.55 ± 0.02 


LfAR,27 


0.55 ± 0.08 


0.56 ± 0.24 


0.62 ± 0.03 



Note. — Lall of Data-MC calculation is performed 
with correcting the OSC luminosity by the overestimate 
of kinematic distance as mentioned in Section 12.2.11 
This correction decreases the luminosity by a ~ 12%. 
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